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We present first measurements of the evolution of the differential transverse momentum correlation 
function, C, with collision centrality in Au+Au interactions at ^/snn = 200 GeV. This observable 
exhibits a strong dependence on collision centrality that is qualitatively similar to that of number 
correlations previously reported. We use the observed longitudinal broadening of the near-side peak 
of C with increasing centrality to estimate the ratio of the shear viscosity to entropy density, n/s, 
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of the matter formed in central Au+Au interactions. We obtain an upper limit estimate of tj/s that 
suggests that the produced medium has a small viscosity per unit entropy. 



PACS numbers: 25.75.Gz, 25.75.Ld, 24.60. Ky, 24.60.-k 

Measurements carried out at the Relativistic Heavy Ion 54 
Collider (RHIC) during the last decade indicate that a 55 
strongly interacting quark gluon plasma (sQGP) is pro- 56 
duced in heavy nuclei collisions at very high beam en- 57 
ergies 1\ . It has emerged that this matter behaves as a 58 
"nearly perfect liquid" , i.e., a fluid which has a very small 59 
shear viscosity per unit of entropy [U [2] • It is a fascinat- 
ing observation that the medium produced in relativistic 
heavy ion collisions reaches exceedingly large tempera- 
tures, of the order of 2 x 10 12 K [3i , in stark contrast to 60 
the very low temperature, T < 3 K, required to achieve 
superfluid 4 He [4]. 61 

Conclusions concerning the shear viscosity per unit of 62 
entropy of the medium produced in Au+Au collisions at 63 
RHIC are based largely on comparisons of non-dissipative 64 
hydrodynamical calculations of the time evolution of col- 65 
lision systems with measurements of the particle produc- 66 
tion azimuthal anisotropy characterized by the elliptic 67 
flow coefficient v-i [H [5] . These calculations describe the w 
V2 and momentum spectra measured in Au+Au collisions 69 
at ytspfN — 200 GeV well at midrapidity (|?7| < 1.0), low 70 
transverse momentum (px < 1 GeV/c), and for mid- 71 
central collisions (impact parameter b < 5 fm) [TJ [6] . 72 
A measure of fluidity is provided by the ratio of shear 73 
viscosity, 77, to entropy density, s, henceforth referred to 74 
as rj/s. It has been conjectured that the limit for all 75 
relativistic quantum field theories at finite temperature 76 
and zero chemical potential is close to the Kovtun-Son- 77 
Starinets (KSS) bound, r}/s\ KSS = (4tt) _1 « 0.08 [3 [7]. 78 
Estimates of r/fs based on V2, measured in Au+Au col- 79 
lisions at ^/snn — 200 GeV, range significantly below so 
the viscosity per unit of entropy ratio of superfluid 4 He ai 
and very close to the quantum limit [2, 5, 8, 9J. Given the 82 
importance of viscosity in furthering our understanding 83 
of QCD matter, it is of interest to consider alternative 84 
measurement techniques to estimate the magnitude of as 
7]/s. Measurements of di-hadron correlations in heavy 86 
ion collisions, carried out as a function of the relative 87 
azimuthal particle emission angle, A</>, have greatly ad- 88 
vanced the studies of hot and strongly interacting matter 89 
at RHIC 10 j. Indeed, studies of correlations between low 90 
and high px particles have revealed the modification of 91 
away-side (A<fi ~ ir) jets and the formation of a longitu- 92 
dinally elongated near-side (A(f> ~ 0) structure, known as 93 
the ridge, in central Au+Au collisions [UJ. Meanwhile, 94 
low-px di-hadron correlation studies reveal rich correla- 95 
tion structures, particularly on the away-side [TT]. How- 96 
ever, the interpretation of these different measurements 97 
is nontrivial, and a number of competing models invok- 98 
ing different reaction mechanisms have been suggested to 99 
explain the data, each with relative success [13J . Thus, 100 



additional observables and measurements are required to 
discriminate fully among these competing models. 

In this work, we present measurements of the differen- 
tial extension of an integral observable C [5] in Au+Au 
collisions at ^Jsnn = 200 GeV. The correlation function 
C is defined as follows: 

/ "1 "2 \ 

(EE PT,iPT,j) 

C(A v ,A<j ) ) = V =i y =1 _X_ (pr)i(pr ) 2 (i) 

where (px)k = (EPT,i) fc / (n) k 1S tne average momen- 
tum, the label k stands for particles from each event and 
the brackets represent event ensemble averages. (n) fe 
is the average number of particles emitted at (rjk,4>k)- 
The indices i and j span all particles in a (rjk,(j>k) bin. 
Ar/ = 771 — r)2 and A<f> = 4>\ — 4>2 are the relative pseudo- 
rapidity and azimuthal angle of measured particle pairs, 
respectively. 

The correlation observable C(Arj, Acp), defined above, 
is an extension of the number correlation function R2 
used in various studies 14 1. By construction, it mea- 
sures the degree of correlation between particles emitted 
at fixed relative pseudorapidity, A77, and azimuthal an- 
gle difference, A<p, and is as such sensitive to various 
aspects of the collision dynamics. However, the explicit 
transverse-momentum weighing provides for additional 
sensitivity to discriminate and study soft (low pt) vs. 
hard (highp^) processes. Note that C differs structurally 
and quantitatively from the observables (SptSpt) |15j 
and AcTp T [TOj previously reported by STAR. Differences 
stem from the fact that C is sensitive not only to num- 
ber density fluctuations, but also to px fluctuations, and 
as such reflects the magnitude of in-medium momentum 
current correlations [5]. 

This study is based on an analysis of 8 x 10 6 mini- 
mum bias (MB) trigger events recorded by the STAR 
experiment in the year 2004 (RHIC Run IV). The MB 
trigger was defined by requiring a coincidence signal of 
two zero-degree calorimeters (ZDCs) located at ±18 m 
from the center of the STAR Time Projection Cham- 
ber (TPC). Data were acquired with forward (+z— axis) 
and reverse {—z— axis) solenoidal magnetic field polarity 
with nominal field strength of 0.5 T. Collision centrality 
was estimated based on the uncorrected primary track 
multiplicity within \r)\ < 1.0. Nine centrality classes cor- 
responding to 0-5% (most central), 5-10% up to 70-80% 
(most peripheral) of the total cross-section were used. A 
mean number of participants, N part , is attributed to each 
fraction of the total cross-section using a Glauber Monte 
Carlo simulation [TT] . 
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1 The analysis is restricted to charged-particle tracks 

2 measured in the TPC with \r/\ < 1.0. Particles of interest 

3 for our measurement are those emerging from the bulk of 

4 the matter. Comparisons of RHIC data to hydrodynamic 

5 models show that the (near) equilibrium description only 

6 holds for particles with pr < 2 GeV/c. For larger mo- 

7 menta, particle production is dominated by hard pro- 
s cesses. Thus, we restrict this measurement to low px, i.e., 

9 with both particles in the range 0.2 < p T < 2.0 GeV/c. 

10 Tracks were selected on the basis of standard STAR qual- 

11 ity cuts |18j . To minimize acceptance effects, events were 

12 analyzed provided their collision vertex lay within a dis- 

13 tance of \z\ < 25 cm from the center of the TPC. How- 

14 ever, the particle acceptance exhibits a small dependence 

15 on the collision vertex position, which may introduce ar- 

16 tificial correlations in the measurement of C . To avoid 

17 such effects, we measure C independently for forward and 
is reverse magnetic field settings in 20 vertex- z bins of width 

19 Az = 2.5 cm in the range —25 < z < 25 cm. Then we av- 

20 erage these measurements to obtain the correlation func- 

21 tion. Track reconstruction inefficiencies for pairs with 

22 Arj ~ 0, due to track crossing or merging in the TPC, are 

23 corrected for by performing a px and charge sign ordered 

24 analysis of these pairs. Track pair losses occur when two 

25 tracks pass nearby one another and produce overlapping 

26 charge clusters in the TPC. For instance, with a forward 

27 magnetic field setting (i.e. along the +z-axis), two posi- 

28 tive charged particles, with pt.i > PT.ii may cross in the 

29 TPC if emitted at pseudorapidity difference A77 ~ 0, and 

30 relative angle A(j> < thereby resulting in pair losses for 

31 A(f> < 0. Pairs emitted with A(j> > however tend to di- 

32 verge in the TPC and thus are not subject to such losses. 

33 In symmetric A+A collisions, pair correlation functions 

34 are invariant under A<j) — ¥ —A(f> reflection. The lost pair 

35 yield at A0 < may thus be corrected based on the 

36 yield at — A<j). Same-sign track pairs are recorded with 

37 Acf> = — \A<j>\ for pt,i > Pt,2 and A</> = +|A0| otherwise. 

38 Pair yields measured for — 1.0 < Acf> < 0, are then substi- 

39 tuted for those at < A(f> < 1.0, thereby compensating 

40 for pair losses. A similar technique is used for unlike-sign 

41 pairs. However, no corrections are made for track pairs 

42 with \Arj\ < 0.032 and \A<f>\ < 0.087 radian (bin at the 

43 origin) . These corrections change the amplitude of C by 

44 < 1% in peripheral collisions and up to 4% in central 

45 collisions. The measurements of C(Ai], A(f>) reported in 

46 this work were constructed using 31 and 36 bins along 57 

47 the A77 and Acf> axes respectively. We verified that the 58 

48 results are independent of the bin width. 59 

49 Figure [T] presents the correlation function, C, for three 60 

50 representative collision centralities (a) 70-80%, (b) 30- ei 

51 40% and (c) 0-5%. Relative statistical errors range from 62 

52 0.8% in peripheral collisions to 0.9% in the most cen- 63 

53 tral collisions at the peak of the distribution. Sources 64 

54 of systematic errors on the amplitude and shape of the 65 

55 correlation function include the collision centrality defi- 66 

56 nition on the basis of primary particle multiplicity in the 67 



(a) 70-80% 




FIG. 1: (Color online) Correlation function, C, shown for (a) 
70-80%, (b) 30-40%, and (c) 0-5% centrality in Au+Au colli- 
sions at -^Siviv = 200 GeV. Cis plotted in units of (GeV/c) 2 , 
and the relative azimuthal angle A<f> in radians. 



range |?7| < 1.0, finite centrality bin width effects, loss of 
track reconstruction efficiency at px < 0.5 GeV/ c, B-field 
direction, and high TPC occupancy, as well as contami- 
nation of the correlation function from weakly decaying 
hadrons (Kg, A), conversion electrons, and HBT corre- 
lations. A study of the effect of the centrality defini- 
tion based on particle multiplicity in the range \rj\ < 0.5, 
|?7| < 0.75, and < 1.0 compared to that obtained with 
the ZDC energy reveals that the \r]\ < 1.0 based centrality 
definition least biases the shape of C at large A-q. Uncer- 
tainties on the correlation yield associated with centrality 
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FIG. 2: (a) Projection of the correlation function C, for 
\A(j)\ < 1.0 radians on the Arj axis for 70-80% centrality, (b) 
30-40% centrality, and (c) 0-5% centrality in Au+Au collisions 
at ysjvjv = 200 GeV. The correlation function C is plotted 2 8 
in units of (GeV/c) 2 . The solid line shows the fit obtained 29 
with Eq. [2] The dotted line corresponds to the baseline, 30 
b, obtained in the fit and shaded band shows uncertainty in 
determining b. 



boundaries and bin width vary from 10% in peripheral 
to less than 1% in the most central collisions. Contam- 
ination from weakly decaying particles and conversion 
electrons is estimated to contribute less than 2% based 
on measured yields and known material budget of the 
detector. HBT effects are essentially negligible, due to 
the large pr range used in the measurement. 

The overall strength of C decreases monotonically from 
peripheral to central collisions. In 70-80% peripheral col- 
lisions, C exhibits a near-side peak centered at A<f> ~ 
A77 ~ and a longitudinally extended away-side struc- 



ture (i.e., broad in Ary) at A</> ~ it. This away-side struc- 
ture largely results from effects associated with momen- 
tum conservation [19]. In more central collisions, momen- 
tum conservation effects are diluted by increased particle 
multiplicities, and the near- and away-side observed cor- 
relation features may result from a superposition of sev- 
eral mechanisms possibly including resonance and cluster 
decays, radial flow effects, anisotropic flow effects, ini- 
tial state fluctuations, and modified jet fragmentation. 
In mid-central collisions (30-40%), the correlation func- 
tion exhibits a sizable broadening of the near-side peak 
and the formation of a near-side ridge-like structure, as 
well as a strong elliptic flow, cos(2A0), modulation [20j . 
In the most central collisions (0-5%), we observe further 
longitudinal broadening of the near-side peak while the 
cos(2A0) modulation and away-side structures have a 
much reduced amplitude. 
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FIG. 3: RMS as function of the number of participating nucle- 
ons for the correlation function C, for nine centrality classes 
in Au+Au collisions at ^/sjvjv = 200 GeV. The dotted line 
represents a lower limit estimate of the RMS explained in the 
text and the shaded band represents systematic uncertainties 
on the RMS. 

We next focus on the longitudinal broadening of C 
with increasing N part based on A77 projections in the 
range |A0| < 1.0 radians. Figures [2ja-c) show the pro- 
jections for 70-80%, 30-40%, and 0-5% centralities, re- 
spectively. The dip seen at Arj ~ for 0-5% central 
collisions (Fig. |2jc) ) is a consequence of track merging 
occurring at A</> ~ Aij ~ 0. We observe that the shape 
and particularly the width of the projections evolve with 
collision centrality. We characterize the widths of the dis- 
tributions by calculating their RMS above a long range 
baseline, b, assumed to be constant in the acceptance of 
our measurement. The baseline, b, is determined using 
the following ansatz to fit the projections: 



g(b,a w ,(j w ,a n ,cr n ) = b + a w exp (- Ar? / 2al, 
+a n exp (-A?7 2 /2cr2) 



(2) 



43 where a w and a n stand for the amplitude of wide and 
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1 narrow Gaussians with widths <r w and cr„, respectively. 57 

2 The offset, narrow Gaussian, wide Gaussian, and full fit 58 

3 are shown in Fig. [2ja) for peripheral collisions. The fits 59 

4 have \ 2 P er degree of freedom values of order unity. The eo 

5 fits are used uniquely for the determination of the offset 61 

6 b. The amplitudes and widths of the Gaussians are not 62 

7 used in the remainder of this analysis. Uncertainties in 63 

8 the determination of the offset, 6, are shown as dark gray 64 

9 shaded areas in Fig. [2j 65 

10 Figure [3] shows the RMS of the correlation function 66 

11 as a function of N part . Vertical lines indicate statisti- 67 

12 cal errors whereas systematic uncertainties on the RMS 68 

13 are indicated by the gray shaded band. Systematic un- 69 

14 certainties arise from several sources. The correlation 70 

15 width exhibits small instrumental dependencies on the 71 

16 magnetic field direction, and the collision vertex position 72 

17 of the order of 3% and 4% respectively in most central 73 
is collision and much smaller in peripheral collisions. Track 74 

19 merging corrections, discussed above, account for parti- 75 

20 cles losses at \Arj\ ~ 0, \A<j>\ < 1.0 and lead to negligible, 76 

21 -C 1%, systematic errors on the RMS of the distribu- 77 

22 tions. The correction technique used does not account 78 

23 for losses at \Arj\ < 0.032 and \A<p\ < 0.087 radian (bin™ 

24 at the origin) which are most severe in 0-5% central col- 80 

25 lisions. This bin is also subject to contamination from 81 

26 e + e~ pairs resulting from photon conversions within the 82 

27 apparatus. We estimated the latter two effects introduce 

28 small systematic uncertainties, < 2%, on the RMS of the 83 

29 correlation functions. The largest source of systematic 

30 uncertainties stems from the baseline determination and 84 

31 the lack of knowledge of the correlations long A77 range 85 

32 behavior, particularly in central collisions. In order to 86 

33 study these effects, we first estimated a lower bound of 87 

34 RMS values, shown as a dotted line in Fig. 3, by set- sa 

35 ting the offset equal to the value of the correlation sig- 89 

36 nal at A 77 = 2.0. This simplistic calculation shows that 90 

37 the RMS exhibits a monotonic growth from peripheral to 91 

38 central collisions. In peripheral collisions, the correlation 92 

39 peak stands atop an approximately flat background but 93 

40 in most central collisions the peak is manifestly broader 94 

41 than the acceptance and this simple estimate is therefore 95 

42 incorrect. We thus used Eq. [2] and systematically studied 96 

43 fits for various number of parameters and fit ranges. Es- 97 

44 timated systematic uncertainties on the offset are shown 98 

45 as gray bands in Fig. 2. Uncertainties on the offset and 99 

46 shape of the distribution, particularly in central colli-100 

47 sions, lead to systematic uncertainties on the RMS rang-101 

48 ing from 10 % in peripheral collisions to 15 % in mostio2 

49 central collisions. The above systematic uncertainties areio3 

50 summed in quadrature and shown as a gray shaded band™ 

51 in Fig. 3. The RMS exhibits a modest increase in the™ 

52 range N part < 100 which may in part result from long™ 

53 range multiplicity fluctuations and from incomplete sys-™ 

54 tern thermalization achieved in small collision systems. ™ 

55 The RMS rises rapidly in the range 100 < N part < 250™ 

56 after which it levels off. no 



According to [8 , shear viscosity should dominate the 
broadening of the correlation function for sufficiently 
large and nearly thermalized collision systems. It should 
thus be possible to utilize the observed broadening to 
estimate the viscosity of the matter produced in these 
collisions. However, jets and jet quenching could also in 
principle contribute to changes in the shape and broad- 
ening of the width of the correlation function with vary- 
ing collision centralities. To examine this possibility, we 
repeated our analysis in the 0.2 < px < 1.0 GeV/c 
and 0.2 < pt < 20.0 GeV/c ranges. Our study shows 
that particles accepted between 0.2 < pt < 20.0 GeV/c 
produce essentially identical widths in peripheral colli- 
sions. In central collisions, RMS reduces by ~7% from 
the RMS widths obtained for the px selection 0.2 < p T < 
2.0 GeV/c. However, lowering the upper pt cut to 1.0 
GeV/c (0.2 < p T < 1.0 GeV/c) does not change the 
widths within statistical errors for 0.2 < pt < 2.0 GeV/ c 
range for the most central collisions, and decreases the 
widths by ~10% in peripheral collisions. We conclude 
that broadening effects associated with jets or jet quench- 
ing are thus likely limited to less than a 10% effect on the 
RMS from peripheral to central collisions. We thus pro- 
ceed to estimate the shear viscosity per entropy of the 
matter produced in central collisions based on the fol- 
lowing formula from Ref . [5] . 

■;-S-^(f 1 -'d) w 

where a c and ag stand for the longitudinal widths of the 
correlation function in central collisions and at forma- 
tion time, respectively, tq refers to the formation time 
and t c t is the kinetic freeze-out time at which particles 
have no further interactions [5T]. T c stands for a charac- 
teristic temperature of the system through its evolution, 
and is here taken to be the critical temperature. We pro- 
ceed by assuming that viscous broadening dominates the 
increase in C with increasing centrality observed in this 
analysis and utilize Eq. [3] to estimate 77/s. We estimate 
cr = 0.54 ± 0.02(stat.) ± 0.06(sys.) by extrapolating the 
RMS width of C to N part ~ 2. The RMS value for most 
central collisions is a c = 0.94 ± 0.06(stat.) ± 0. 17 (sys.). 
Using commonly accepted estimates of 1 fm/c, 20 fm/c, 
and 170 MeV [5j for the formation time, central collision 
freeze-out, and effective temperature, we obtain a value 
of r)/s = 0.13 ± 0.03. Inclusion of systematic uncertain- 
ties on the widths leads to a range of 77/s = 0.06 — 0.21. 
Figure |4] shows 77/s as a function of Tq 1 — t^ 1 and pro- 
vides an estimate of theoretical uncertainties based on a 
literature survey of theoretical estimates for To and T c . tq 
is typically assumed to be in the range 0.6 - 1.0 fm/ c (e.g., 
[H |2lJ [22]). Here, we have assumed that the broadening 
of C is entirely due to viscous effects. Given that other 
(unknown) dynamical effects could perhaps also lead to 
the correlation function broadening, we conclude that our 
measurement provides an upper limit. Based on the sta- 
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tistical and systematic uncertainties of our measurement 37 
(using one standard deviation) and caveats of the used 38 
theoretical model, and using the ranges 150 < T c < 190 39 
MeV and 0.6 < t^ 1 — < 1.6 (fm/c) -1 , we derive an « 
upper limit of order rj/s ~ 0.3 . « 




xJ-t-i [(fm/c)" 1 ] 



FIG. 4: 77/s as a function of r 1 — r c j and T c . tq and T c 
vary from 0.5 < r < 1.5 fm/c and 150 < T c < 190 MeV, [ 
respectively. t 

In summary, we presented first measurements of the e 
differential transverse momentum correlation function C t 
from Au+Au collisions at ^Jsmn — 200 GeV. In pe- 
ripheral collisions, C has a shape qualitatively similar to t 
that observed in measurements of number density cor- t 
relations, with a relatively narrow near-side peak near e 
Arj « A(/3 « and a longitudinally broad away-side e 
[Ttil US] . We find that the near-side peak progressively e 
broadens with increasing number of collision participants e 
while the overall strength of the correlation function de- _ 
creases monotonically. These results may be used to fur- 7 
ther constrain particle production and correlation mod- ? 
els. We used the observed longitudinal broadening to es- 
timate 77/s of the matter formed in central Au+Au colli- 7 
sions. Considering systematic uncertainties in the deter- ' 
ruination of correlation widths, particularly in central col- ' 
lisions, and assuming somewhat conservative estimates _ 
of the temperature, formation and freeze-out times, we E 
obtain a range of 77/s = 0.06 - 0.21. This result is remark- £ 
ably close to the KSS bound, (47r)~ , and is consistent £ 
with results obtained from hydrodynamical model com- £ 
parisons to elliptic flow data [5]. 
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